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ABSTRACT 



Using the longest optical-interferometeric baselines currently available, we have 
detected strong near-infrared (NIR) emission from inside the dust-destruction 
radius of Herbig Ae stars MWC275 and AB Aur. Our sub-milli-arcsecond resolu- 
tion observations unambiguously place the emission between the dust-destruction 
radius and the magnetospheric co-rotation radius. We argue that this new com- 
ponent corresponds to hot gas inside the dust-sublimation radius, confirming 
recent claims based on spectrally-resolved interferometry and dust evaporation 
front modeling. 

Subject headings: young stellar objects — circumstellar disks — radiative transfer 
— dust sublimation — interferometry 



1. Introduction 



The inner tenths of an AU in disks around young stellar objects (YSOs) hosts fascinating 
astrophysical phenomena like disk winds, jets, dust sublimation/condensation and thermal 
annealing, and gas in-fall onto the star. To make spatially resolved observations of these 
phenomena, milli-arcsecond angular resolution is needed even f or nearest the systems. Using 
the sub-milli-arcsecond resolution capabilities of the CHARA (jten Brummelaar et al.ll2005l ) 
interferometer array (operated by Georgia State University on Mt. Wilson), we probed inner 
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disks of the YSOs MWC275 and AB Aur at 0.1AU scales and found NIR emission within 
the dust destruction radius in these systems. 

MWC275 (Ale, d=122pc, V=6.86, K=4.59, M*=2.3M ) and AB Aur (AOpe, d=144pc, 
V=7.01, K=4.37, M*=2.4M Q ) are Herbig Ae stars - pre-main-sequence stars of i ntermediate 
mass . This class of stars shows a significant excess of near-infrared (NIR) flux (INatta et al. 



2001 ) over the stel l ar photospheric emission . Current theories (iDullemond et al 



2001 



Isella fc Nattal 120051 ; iTannirkulam et al.l 120071 1 explain this excess solely on the basis of a 
"puffed up inner dust rim". In these "rim" theories, dust sublimation terminates the cir- 
cumstellar dust disk at a finite radius (the truncation radius depends on the luminosity of 



the central star and dus t properties) creat ing an optically-thin cavity (jMillan-Gabet et al. 



19991 ; iTuthill et al.l l200ll ; Eisner et al.l 120041 ) around the star. The outer edge of the cavity 
forms a dusty rim that traps stellar photons, re-radiating the energy mostly in the NIR. 
The rim idea seems to broadly explain much of the correlation b etween the NIR sizes and 



luminosities of the stars observed in a number of Y SO systems (IMonnier fc Millan-Gabet 
20021 : IMonnier et aDbood : iMillan-Gabet et alil20~07V 



However, recent results have begun to indicate that models in which the entire NIR 
excess aris es solely from dust r ims are inadequate in explaining some of the observed YSO 
properties. lAkeson et al.l (120051 ) noted that NIR emission from the the dust rim al one was not 



suffici ent to explain the spectral energy distribution (SED) of some TTauri stars. Ilsella et al. 



(120061 ) pointed out that th e dust evaporation temperatures of well studied astronomical dust 



Eisner et al. 


( 


2007 


) and 


Eisner 


( 


2007) 



had to invoke gas emission to account for the apparent decrease in NIR si zes of Herbig 



Ae/Be stars with decreasing wavelength. Most recently, iKraus et al.l (120071 ) showed that 
the NIR visibilities of Herbig Be star MWC147 could not be explained without including 
emission from inside the dust evaporation radius. Due to the lack of sufficient spatial resolu- 
tion, some of these conclusions are tentative and based on indirect arguments with multiple 
interpretations. 



2. Observations 

To unambiguously constrain NIR disk structure in YSOs, we began an observation cam- 
paign with the CHARA Array using the "CHARA Classic" beam combiner (Iten Brummelaar et al 
20051 ) . The targets were observed at K band (central wavelength of 2.13 microns) with the 
long baselines of CHARA at a variety of orientations in 8 runs between June 2004 and 
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June 2007. The longest baseline observation for MWC275 was 325m (resolution of 0.67 
milli-arcseconds) and 320m (resolution of 0.68 milli- arcseconds) for AB Aur. The data were 
reduced using standard CHARA reduction software (Iten Brummelaar et al.ll2005l ) and these 
results were cross checked with an independent code developed at University of Michigan 
(the full calibrated OIFI TS data set will be p osted on OLBIN - jhttp: //olbin.jpl.nasa.gov] ). 
HD164031 , HD166295 jMerand et all 12005^ and HD156365 with uniform-disk (UD) di- 
ameters of 0.83± 0.08 mas, 1.274±0.018 mas and 0.44±0.06 mas were used as calibrators 
for MWC275. AB Aur visibilities were calibrated with HD29645 (UD diameter=0.54±0.07 
mas) and HD31233 (UD diameter =0.76±0. 13 mas). The visibility errors (which include 
calibration errors, statistics and errors due to finite calibrator size) are at the ~6% level, 
typical for CHARA Classic. The ne w fringe-visibility d ata fro m CHARA were c ombined 
with past measureme nts from IOTA (IMonnier et al.ll2006l ). PTI (jEisner et al.l 120041 ) and the 
Keck Interferometer (IMonnier et al.l 120051 ) providing a baseline coverage from 15m-320m on 
the targets, placing strong constraints on the geometry of the inner-most disk regions. 



The baseline orientations used have allowed us to clearly detect the asymmetry of the 
MWC275 disk (details in a followup paper), as having inclination=48°±2°, P A=136°±2°, 
consis tent with the inclination of Sllg 1 degrees, and PA of 139°±15° determined in lWassell et al 
(120061 ). Thus, in Fig $2]the visibility of MWC275 is shown as a function of "effective baseline" 



B. 



eff 



Bprojected V COS 2 (9) + COS 2 (0) Sffl 2 (9) 



where 6 is the angle between the uv vector for the observation and the major axis of the disk 
and (p is the inclination of the disk. Under the flat disk assumption, the effective baseline 
correctly accounts for the change in resolution due to the disk inclination and PA (geometry 
of finitely thick disks is represented only approximately with optical depth effects and 3- 
D geometry of thick disks not being taken into account), allowing us to plot the visibility 
measurements as a function o f one coordinate, s implifying presentation and analysis. As the 
AB Aur disk is near face on (jEisner et al.l 120041 ). we have plotted its visibility as a function 
of true projected baseline (Fig §2]). 



L Resolution is denned as ^ -lambda is wavelength of observation and D is the interferometer baseline 



length. 
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3. Modeling 



To elucidate the structure of the inner disk and the dust rim, we constructed dust-only 
rim models for MWC275 and AB Aur and fit the models to the visibility data. In these 
models, all of the K-band emission is assumed to come from the inner dust rim and the 
central star ("standard" model). The rim is in hydrostatic equilibriu m and its shape is set 



by the pressure dependen ce of the dust grain evaporation temperatures (jlsella fc Nattall2005 



Tannirkulam et al.ll2007l ). The normalization of the pressure dependence of the evaporation 
temperatures (discussed in more detail in a followup paper) is adjusted so that the model 
rims fit the short baseline (< 100m) visibility data. As can be seen in Fig §TJ a realistic 
treatment of dust physics produces a curved inner rim with a sharp inner edge, set by dust 



destr uction, and a sharp outer edge, set by the fact that the rim shadows (IDullemond et al. 



20011 ) a portion of the disk behind it preventing direct star-light from reaching the region. 



We are able to fit the standard model to the visibility data for MWC275 and AB Aur 
at baselines shor ter than 100m (Fig jj2D fai r ly wel l [ f| (a result already e s tablished in multi ple 
previous studies; iMonnier fc Millan-Gabetl J2002h : feisner et all J2004h : llsella et all J2006h ). 
Rim radii from the model fit are 0.19 AU (thickness 0.04 AU) for MWC275 and 0.21AU 
(thickness 0.05 AU) for AB Aur (see Table £Q] for details). However, the standard model 
fails glaringly at fitting the long-baseline data. For baselines longer than 110m the standard 
rim predicts strong bounces in the visibility response which are not observed in the data. 

The lack of bounce in visibility cannot be explained in the framework of the standard 
rim model. Models for NIR excess in Herbig Ae stars in which all of the excess arises in rims 
will display bounces in visibility at long baselines because of the presence of sharp ring-like 
features with high spatial frequency components in the corresponding images, even for the 
smoothest rims physically plausible. 

The only explanation for the absence of long-baseline visibility bounce in the data is 
the presence of an additional smooth emission component within the dust destruction radius 
that suppresses the large visibility bounces seen in a rim-only model (the smooth emission 
component has to be on size scales similar to the dust sublimation radius and cannot arise 
from a compact source like magnetospheric accretion). Fig fjl] shows a standard model to 
which NIR emission has been added within the dust evaporation radius. Our data cannot 



2 The near-IR disk size obtained from the Keck Interferometer data is ~20% larger than the size obtained 
with the CHARA data (June2004-Aug2006 epochs). The size determined from the S2W1 June2007 data also 
differs at the ~25% level from the size obtained from earlier CHARA epoch s. The size fluctu ations could 
be related to the observed variability in the infrared photometry of MWC275 (jSitko et al.ll2007l ) . Variability 
does not affect any of our conclusions. 
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uniquely constrain this emission profile and we assume a constant surface brightness disk (a 
uniform disk), which is a simple model. A rim + uniform disk model in which the uniform 
disk emission (refer Table £JT] f° r details) accounts for 40-60% of the total K-band emission 
provides a good fit to the visibility (Fig §2]) over the whole baseline range. 

Besides the interferometry, emission inside the dust evaporation radius can also help 
explain the NIR flux budget of MWC275 and AB Aur. MWC275 and AB Aur have total 
K-band emission a factor of 10 larger than the stellar p hotospheric K-band e mission, whereas 



the total emission produced by standard rim models (jlsella fc Nattall2005l ) is only a factor 
of 3.5 larger (SED analysis will be presented in a followup paper) than the photospheric 
emission. The emission component inside the dust rim can supply the "missing" K-band 
flux in MWC275 and AB Aur. In the next section we discuss candidates for the physical 
origin of this emission. 



4. Physical Interpretation 



We consider 3 opacity candidates for pr oviding the smooth K -band emission component: 
a) a dusty spherical halo around the stars (jVinkovic et al.ll2006l ). b) highly optically trans- 
parent /refractory dus t extending inside the rim and c) gas inside the dust evaporation front. 
Vinkovic et al.l (120061 ) showed that halo models can fit NIR SED and some early NIR inter- 
ferometry data on Herbig Ae stars. However, the detecti on of asymmetries in NIR emission 
of MWC275 (Section §2]) and other Herbig Ae/Be stars (lEisner et al. ||2004h support a disk 



model over the spherically symmetric halo model of IVinkovic et al 



(120061) . The gas densitie s 



in the inner disk are not high enough to let dust survive beyond 1800K (jPollack et al.lll994l ). 
making optically transparent /refractory dust an unlikely candidate for the excess K-band 
emission. 

Hot gas is an excellent candidate as a NIR emission source inside the dust sublimation 
radius . Molecular gas is known to have den se NIR opacity features at kilo-kelvin tempera- 
tures (ICarr et al.ll2004l ; iFerguson et al.ll2005l ). The required gas emission levels to explain the 
long-baseline visibility data can be obtained with optic ally thin gas emission (optical dept h 
t ~ 0.15) with a temperature range of 2500K-3500K jMuzerolle et alil2004! ; Eisneri bpOTt ). 
Our claim regarding the presence of gas emission inside t he dust rim is furt her supported 
by the detection of the Br7 line in young star HD 104327 (ITatulli et al.l 120071 ) on size scales 
similar to the dust destruction radius. 



The required color excess (E(H-K)) from gas emission to match the SED in MWC275 
and AB Aur is ~ 0.1. MWC275 and AB Aur have Br7 emission equivalent widths of 
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~6A ( Garcia Lopez et al.l 12006 ). The equivalent widths are in principle sufficient to explain 
the color excess with free-free emission based on the observed correlation between color 



excess and Br7 equivalent width for Be stars (IHowells et al.ll200ll ). Thus free-free emission 
is also a good candidate for NIR emission inside the dust evaporation front. 

At first thought, presence of binary companions or source variability would also seem 
good explanations for the excess K-band emission and the lack of visibility bounce at long 
baselines. However, the presence of a binary companion within 0.5 AU of MWC275 that 
contributes ~50% of the total near-IR emission will cause changes of ~0.5 in the absolute 
visibility of the system at any given baseline over the multiple CHAR A observation e pochs . 

320061 ) 



Furthermore, Monnier et al 



This variation is not seen in the data (see Fig. §! 
obtained NIR visibility and closure-phase data on MWC275 and found no evidence for a 
binary companion ruling out binarity as an explanation for the observations. A comparison 
between the Keck Interferometer and CHARA data (taken in different observation epochs) 
in Fig. $2] shows that MWC275 visibilities are time variable at the ~20 % level, suggesting 
that variability could complicate our analysis. However, the data from the S2W1 CHARA 
telescope pair (Fig. ||2]) were obtained on a single night and show a similar visibility-baseline 
trend as the rest of the data, disfavoring variability as an explanation for the observed trends. 
Thus, our conclusion on the presence of smooth gas emission inside the dust destruction 
radius is robust. 



5. Discussion 

In this work, we have conclusively shown the presence of strong NIR emission between 
the magnetospheric co-rotation radius and the dust sublimation radius in Herbig Ae stars 
MWC275 and AB Aur. The underlying mechanism for this emission is not known requiring 
high resolution NIR spectroscopy and self-consistent modeling of the gas-dust transition 
region to identify the gas species and the physical processes responsible for the emission. 

The detection of continuum gas emission within the dust destruction radius through 
long-baseline interferometry is of special importance to understanding the star-disk connec- 
tion. Accretion models for YSOs have tended to focus on two separate physical scales (i) 
the process of mass and angular momentum transport in the outer disk and (ii) gas ac- 
cretion from a few stellar radii on to the stellar photosphere. Milli- meter interferometry 
and direct imaging in the mid-IR have allowed observational access to the outer disk while 
high-resolution optical-spectroscopy can probe the circumstellar environment around YSOs 
on scales of a few stellar radii. Our near-IR interferometry results on MWC275 and AB Aur 
provides exciting prospects for connecting the two scales and determining the structure of 
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the gas flow from a few tenths of an AU down to the stellar surface. In the following para- 
graph, we describe NIR interferometry as a tool in probing accretion physics inside the dust 



rim. 



In the regime where in the inner gas disk is optically thin, the Shakura-Sunyaev viscosity 
parameter for the gas disk can be constrained using a combination of NIR interferometry 
and optical spectroscopy. NIR interferometry in multiple wavelength bands can constrain 
the inner gas-disk surface density and temperature profiles. The surface density of the disk 
i s related to the mass accr etion rate and the viscosity parameter in quasi-steady-state disks 
(jShakura &: Sunyaevlll973l ). Thus, simultaneous measurements of mass in-fall rates (through 
modeling of optical-emission lines and optical con tinuum excess above photospheric emission 
JCalvet fc GullbrineJll998l : iMuzerolle et aDboOlh ) and the gas surface density profiles will 
directly constrain the viscosity parameter 'a' in the inner disk. This is key to understanding 
angular momentum transport in circumstellar disks. 
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Table 1: MWC275 and AB Aur model parameters for K-band emiss ion. The models als o 



include extended envelopes contributing 5% of the K -band emission (IMonnier et al.ll2006f ) 

Star model star dust-rim gas-disk rim 

type emission emission emission inner radius 

(%) (%) (%) (AU) 

MWC275 dust rim only 25 69 - 0.19 

dust rim+gas 10 29 56 0.22 

AB Aur dust rim only 22 73 - 0.21 

dust rim+gas 12 45 38 0.25 
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Standard' inner-rim model (K band image) 




'Standard' inner-rim + gas model (K band image) 




Fig. 1. — Representative models for n ear-IR emission in Herbig Ae stars, a) Left pa nel. A 
standard curved dust-rim-only model (jlsella fc Nattall2005l : iTannirkulam et al.ll2007l) wher e 
the dust sublimation temperature is a power-law function of gas density (IPollack et al.lll994j ). 
b) Right panel. Gas emission (modeled as a uniform disk centered on the star) has been 
added inside the dust rim in order to smooth out the emission profile. 



Visibility curve for MWC275 Visibility curve for AB Aur 




Effective Baseline (m) Baseline (m) 



Fig. 2.— Visibility data and models for MWC275 (left) and AB Aur (right). MWC275 
visibilities are plotted as a function of "effective baseline" (see section $2]), which accounts 
for the change in resolution due to disk inclination and PA. 

The dotted lines correspond to "Standard" rim models tuned to fit visibility data for each 
source for baselines shorter than 100m. The solid lines correspond to dust rim+gas models. 
The model parameters are listed in Table ^TJ 



